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Abstract 

£NJ ■ Motivated by the recent observation of the 511 keV 7-ray line emissions from the galactic bulge 

and an explanation for it by the decays of light dark matter particles, we consider the light axino 
whose mass can be in the 1 — 10 MeV range, particularly, in the context of gauge-mediated su- 
t^. \ persymmetry breaking models. We discuss the production processes and cosmological constraints 

for the light axino dark matter. It is shown that the bilinear R-parity violating terms provide an 

T 1 

' appropriate mixing between the axino and neutrinos so that the light axino decays dominantly to 

e + e~v. We point out that the same bilinear R-parity violations consistently give both the lifetime 
of the axino required to explain the observed 511 keV 7-rays and the observed neutrino masses 
C ! and mixing. 
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I. INTRODUCTION 



The Peccei-Quinn (PQ) mechanism to solve the strong CP problem P, when combined 
with supersymmetry (SUSY) which is the solution to the gauge hierarchy problem, predicts 
a singlet fermion called the axino. It can be light in certain supersymmetry breaking mech- 
anism, and become the lightest supersymmetric particle (LSP) providing a good candidate 
for the particle dark matter (DM) in various mass ranges 0, y, 0, IE El ■ 

Phenomenologically viable supersymmetric models are implemented with the R-parity 
to assure the stability of the proton, which also implies the stability of the LSP. However, 
R-parity is not dictated from any deep theoretical principle. The small violation of R-parity 
is an attractive option for generating the neutrino masses and mixing Even with the 
R-parity violation, the LSP can be cosmologically stable and it may provide an indirect 
detection mechanism of the DM by leaving imprints in 7-rays from the galactic center and 
in the diffuse background 

Recent observation of 511 keV 7-rays by the SPI spectrometer aboard the INTEGRAL 
satellite not only confirmed the previously measured total flux but also revealed the morphol- 
ogy of the bulge emission, which is highly symmetric and centered on the galactic center 
with a full width half maximum of ~ 8° [9|, [lOj, [ll|. The observed emission of 511 keV 



7-rays can be well explained by e + e~ annihilations via positronium formation. But the 
origin of these galactic positrons remains a mystery. Many astrophysical sources have been 
suggested, including massive stars, neutron stars, black holes, supernovae, and X-ray bina- 
ries. The generic problem of astrophysical sources is that they have difficulty in explaining 
both the total flux and the high bulge-to-disk ratio of observed 511 keV 7-rays. Given 
this difficulty, suggested were alternative explanations that light dark matter (LDM) parti- 
cles annihilati ng o r de cay ing in the galactic bulge are the sources of the galactic positrons 

In addition to positrons, annihilations or decays of LDM particles produce direct 7-rays 
via the internal bremsstrahlung processes. The observation of 7-rays from the galactic center 
in the energy range 1 — 100 MeV bounds the mass of LDM particles to be less than about 
20 MeV |l9|. It was also claimed that astrophysical sources are missing for the diffuse 
7-ray background in the energy range 1 — 20 MeV from the observed spectrum, and that 
direct 7-rays from annihilations or decays of LDM particles can fit the spectrum when the 



produced positrons are normalized to fit the 511 keV 7-rays from the galactic bulge |20t 121 
Concerning the annihilating LDM, its mass less than 10 MeV is practically excluded because 
it leads to a much longer supernovae cooling time which makes impossible the emission of 
sufficiently energetic neutrinos observed in SN1987A 2^ . 

In view of above observations, the axino in R-parity violating supersymmetric models is a 
well- motivated candidate for the MeV dark matter whose decay can explain the observed 511 
keV line emission from the galactic bulge as suggested by Hooper and Wang fl3j ]. Indeed, 
R-parity violation is required to make the axino decay and its lifetime can be very long 
since its interactions are suppressed by the PQ scale. An interesting question one may ask 
is whether the same R-parity violation can also generate the observed neutrino masses and 
mixing. 

In this article, we show that the axino LDM scenario is consistent with the usual mecha- 
nism of generating the neutrino masses and mixing at tree-level through the small bilinear 
R-parity violating couplings ~ 10~ 6 j^. Such small bilinear terms turn out to induce an 
appropriate axino-neutrino mixing through which the light axinos decay to the positrons 
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with the right range of the lifetime [1 



4 x 10 26 

Tdm ~ 7U \ a sec • 1 

m dm (MeV) 



This has to be contrasted to the case of [13] where the trilinear couplings Ajii ~ 0.1 were 
considered. 

We also discuss how the MeV axino can arise, particularly, in gauge mediated supersym- 
metry breaking (GMSB) schemes where the saxion is predicted to get the mass in the range 
of 4 — 50 GeV. Axinos are produced thermally or non-thermally in the early universe and the 
amount of axinos can be correctly adjusted for the appropriate reheat temperature and/or 
MSSM parameters. If the saxion abundance is comparable to the axino abundance as is the 
case of the thermal regeneration, the saxion decay to ordinary particles can cause a problem 
of upsetting the standard prediction of the big bang nucleosynthesis (BBN). Such a "saxion 
problem" puts another cosmological constraints on the axino LDM models. 



II. THE AXINO MASS 

The axion supermultiplet A = (s + ia, a) consists of the pseudo-scalar axion a, its scalar 
partner, the saxion s, and its fermionic partner, the axino a. It has the model-independent 
interactions with the gluon supermultiplet W a 

A ff = -^rAWaW* , (2) 

where f a is the PQ symmetry breaking scale. At present particle phenomenology, astrophys- 
ical and cosmological observations restrict the range of f a to be 10 9 GeV ^ f a ^ 10 12 GeV. 
Then the axion mass is given by m a ~ Aq CD // a ~ 10~ 2 — 10 -5 eV. 

The axino mass depends crucially on the way of supersymmetry breaking. In generic 
supergravity (SUGRA) models, it is expected to get the typical soft mass of order m 3 / 2 ~ 100 
GeV and some special arrangement, e.g. no-scale model, is needed to allow the axino mass 
in the MeV scale 0, 0|. Light axino can arise naturally in GMSB models where SUSY 
breaking scale is lower than the PQ symmetry breaking scale 0] • Let us show how the MeV 
axino is predicted in GMSB models. Consider the DFSZ axion model [lj where the MSSM 
fields are charged under the PQ symmetry. Upon the PQ symmetry breaking, an effective 
Kahler potential between the axion supermultiplet A and the other fields $j is generated as 
follows; 



Keff = e** — $1$, (3) 

where Xi is the PQ charge of $j. Taking the terms of order A 2 and $j = H12, one has a 
contribution to the axino mass; m a pa Fh.v/ f 2 ~ [iv 2 ' j ' f 2 <C MeV which is negligible in our 
context. In GMSB models, $j can be one of the hidden sector superfields, say X, which is 
assumed to take the vacuum expectation value; {X) = X + @ 2 Fx leading to the effective 
supersymmetry breaking scale, A = Fx/X = 10 4 — 10 5 GeV |25|]. Then, one obtains the 



axino and saxion mass as 



ma = 4 ^ k " (f) A ' (4) 

m\ = 2x 2 M . (5) 

J a 
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The axino mass in the range 1 - 10 MeV is obtained with X/f a ~ 1(T 3 - 1(T 4 . These 
equations also give us the relation; 

m^2m 5 A (6) 
leading to the saxion mass m s ~ 4.5 — 45 GeV. 



III. THE ORIGIN OF COSMIC AXINOS AND COSMOLOGICAL CON- 
STRAINTS 

There are two known ways in which axinos are produced in the early universe. One is the 
thermal production from the hot thermal bath after reheating. The other is the non-thermal 
production from decays of the lightest ordinary supersymmetric particles (LOSPs). 

The decoupling temperature of axinos is estimated as 

r ^ 101 ° GeV (lo4v)(5T)"- < 7 » 

where a s is the strong coupling constant. If the reheat temperature Tr after inflation is 
higher than the decoupling temperature, the universe is overpopulated by axinos if the 
axino mass is larger than a few keV. Therefore, we only consider the case that the reheat 
temperature is lower than the decoupling temperature. In this case, axinos are produced 
from the thermal bath through scattering of quarks and gluons, though the number density 
of them do not reach the thermal equilibrium. The amount of axinos produced in this way, 
so called regeneration, is estimated to be p 

^ 2s30 ' 28 (^)(lo^v)(lo4v)"- < 8 » 

Thus, for the axino with mass 1 — 10 MeV to be the LDM, the relevant range of reheating 
temperature is 10 — 100 TeV. 

The axinos from decays of LOSPs can be cosmologically interesting when the axino mass 
is around the marginal value of order 10 MeV. For this size of axino mass, the reheat 
temperature must be lower than 10 TeV to suppress the thermal production (regeneration). 
The amount of produced axinos is simply connected to that of LOSPs by 

Sl-tf = H^Oh 2 , (9) 
m x 

and independent of the reheat temperature. When we take m x = 100 GeV and m„ = 10 
MeV, the required value of fi x h 2 is ~ 10 4 . Such a high value is reached for very large MgusY 
in the range of tens of TeV. Thus, the non-thermal production of axinos for LDM could only 
be marginally relevant. 

Even though relic axinos dominantly come from regeneration, the existence of LOSPs 
and their decay to axinos can produce radiative or hadronic cascades during or after the 
BBN, and alter its standard predictions on the light element abundances. To avoid this, 
the mass of LOSP needs to be large enough to make its lifetime much less than 1 sec. For 
example, in the case of the neutralino, one requires m x > 150 GeV. 

Let us discuss here how the accompanied saxion can also upset the standard prediction 
of the BBN, which is called "the saxion problem" j26{. Contrary to the axino, the saxion 
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has the axion-like couplings to the quarks, ^r-sqq, or leptons, ^sll, so that its life-time is 

J a J a 

much shorter than the axino LSP. On the other hand, during the axino regeneration (|8|). 
the saxions are also populated by the same amount and thus one finds 

VMeV/ V m a J V 0.28 / V ; 

where Y s is the saxion number density in unit of the entropy density. Note that this quantity 
is strongly constrained by the BBN. In the mass range m s > O(10) GeV, the above equation 
gives m s Y s > 10~ 5 GeV for mj = 1 MeV. Now that the saxion decays mainly to bottom and 
charm quarks, one finds a strong limit on the saxion lifetime: t s < 10~ 2 sec |22[. Specifically, 
the mass relation © gives us m s w 14 GeV for the axino mass m a ~ 1 MeV and A = 10 5 
GeV. Then, the saxion lifetime, 



1 2 1 - 1 

1 mi 
I — rl m& 



<10" 2 sec, (11) 



becomes short enough to avoid the saxion problem for f a < 3 x 10 11 GeV. In the case of 
supergravity models where one expects to get m s ~ 10 2 ^ 3 GeV, the saxion is free of such a 
problem. 



IV. AXINO NEUTRINO MIXING AND AXINO DECAY 

Let us now assume the generation of the bilinear superpotential term, H1H2, and its 
R-parity and lepton number violating extension, L1H2 as a result of the PQ symmetry 
breaking; 

W eff = nHiH 2 + unUHz (12) 

where \x and 6ifi carry PQ charges whose sizes are determined by the PQ charge assignments 
for Hi 2 and L 4 . In Eq. (j3J), the leading terms in A, 

Keff = Y^HjHi + x L] L)L 3 ] + • • • , (13) 

la 

give rise to the following axino-Higgsino and axino-neutrino mass terms; 

>C m i x i n g = x^—aHx + x H2 — aH 2 + x Li ——aUi + h.c. . (14) 

J a J a J a 

For jiv/ f a <C and e^v^j f a "C m a , one has the approximate mixing angles between the 
axino and Higgsino or neutrino as follows; 

Q~aH ~ J" and 0^ = *Li-f— ■ ^> 

J a J a'' 

The axino-neutrino mixing derived above induces the effective vertex of aviZ and a/jW 7 
with the coupling ~ g9a Ui - This gives rise to the four-quark operator as follows: 

G 

£ e + e - « 0&/< ^7(u75« e7 M (25 a - 7 5 )e (16) 
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where we omitted the small correction due to the vector part of the charged current. 

Another important interaction to consider is the axino-photon-neutrino vertex arising 
from the photino-neutrino mixing. The bilinear term L1H2 induces the mixing between 
neutrinos and neutralinos of order £j. Then the supersymmetric anomaly coupling of axino- 
photon-photino leads to the axino-photon-neutrino coupling which is written down schemat- 
ically as follows; 

= C^Oja ^ (17) 

where C a77 is an order-one parameter taking into account the precise values of the U(l) em 
anomaly and the photino-neutrino mixing. From the vertices (J16|) and (|17|) . we get the 
following decay widths of the axino; 

G 2 F m\ 2 i-l 1 r 1 
Ul£+e ~ ~ 192^3 ^4 ilJ 
C 2 a 2 m 3 

(16vr)3 f a ^- ^ 

Let us first note that the photon mode is suppressed by a 2 m compared to the e + e~ mode; 

r 0^2 2 

Y ve+e - ~ 32G>V ~ [iJ) 

for yu/C a77 = 100 GeV. It is smaller than the internal bremsstrahlung process of e + e~ mode 
which is suppressed by a em and also produces a direct 7-ray. This is enough to be consistent 
with the observations of the MeV 7-ray spectrum [19]. Then, the axino decay is determined 
by the process a — > ue + e~ whose lifetimes is given by 

which is in the right range to explain the observation ([TJ consistently with the neutrino data 
as will be shown in the following section. 



V. CONSISTENCY WITH THE NEUTRINO DATA AND EXPERIMENTAL SIG- 
NATURES 

One of the interesting aspect of R-parity violation is that it can be the origin of the 
observed neutrino masses and mixing [23(. The general superpotential allowing R-parity 
and lepton number violation includes the following bilinear and trilinear terms; 

Wr p = eifiLiH 2 + -Xij k LiLjEl + X' ijk LiQjDl . (21) 



According to the observation of Ref. [Lj|, an appropriate life time of the axino decay a — > 
v ii,r€ + e~ can arise with trilinear R-parity violating couplings A2ii,3n ~ 0.1. Such trilinear 
couplings can generate the 2-3 components of the neutrino mass matrix; 

«5 « ^n^n 1 ^^ (22) 
J 8tt 2 m~ 
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where tan/? = (H®) / (H®) and mg is the selectron soft mass. While the charged- current and 
e-fi-T universality put the bound Xm ^ 0.1(mg/200 GeV) |28!], the above one-loop mass 
can reach the observed atmospheric neutrino mass scale m v pa 0.05 eV only for an extreme 
value of fitan/3 ~ 50 TeV taking the boundary value of Aju = 0.1 (mg/200GeV). In order 
to generate the other components of the neutrino mass matrix, one needs to introduce 
some other trilinear couplings such as A i2 2,j33 which induce M^,M^ 2 and M" 3 through the 
combinations of XijjXijj, A133A233 and Ai 2 2A 32 2, respectively. Then, appropriate neutrino 
masses can be obtained for the trilinear couplings, Aj 2 2 ~ 10 -4 and Aj33 ~ 10~ 5 , where 
the small ratios Xi22/Xm and X^/Xm are dictated by the factors of m e /m M and m e /m T , 
respectively. Such a hierarchy among Xijj appears ad-hoc considering the usual hierarchy in 
the quark and lepton Yukawa couplings. 

Nevertheless, if there exits the trilinear coupling Xm of order 0.1, they leads to a re- 
markable experimental signature of resonant single sneutrino production in the future linear 
collider |29l . l30l | , non-observation of which would rule out the axino LDM decaying through 
the trilinear couplings. 

The observed neutrino masses and mixing can be more naturally explained if one invokes 
the presence of the bilinear term of the order 10~ 6 [23]. The bilinear R-parity violation 
generates neutrino masses at tree-level through the neutrino-neutralino mixing. In addition 
to the 6i term in the superpotential (jUJ), the scalar potential also contains the R-parity 
violating bilinear soft terms as follows; 

V = ml iHl UHl + BiUHi + h.c, (23) 

where Bi is the dimension-two_ soft parameter. Generically, one has Bi = eiBfi with a 
dimension-one soft parameter B for the \x term, and the soft mass-squared m 2 L H contains 
the supersymmetric term e^/i 2 . Upon the electroweak symmetry breaking, the sneutrino field 
gets nontrivial vacuum expectation value; 

(fy _ m 2 LiHi + Bjtanp ^ 



Vi m 



"1 



which is expected to be of order e< up to the soft mass dependence. These bilinear parameters 
induce mixing between neutrinos and neutralinos. For the small mixing mass, the week-scale 
seesaw with heavy neutralino mass scale ~ 100 GeV leads to the well-known neutrino mass 
matrix at tree-level; 

M 2 

M^. = --^^cos 2 /3 (25) 
r N 

where & = a - (i>i)/vi and F N = .U 2 .".U, + M|cos2/3//i with = c^Mi + s 2 v M 2 . 



From Eq. (|25|h one obtains the size of |^| = \£,i\ 2 consistently with the atmospheric 
neutrino mass scale as follows; 

|e| = 0.7xlO-^f^V /2 f^-V /2 . (26) 
M cos/3 \M Z J V0.05eVy V ; 

This is compatible with the axino lifetime relation (J20|) for ^ ~ e^. Note that the smaller 
neutrino mass explaining the solar neutrino oscillation can come from one-loop diagrams 
involving the trilinear couplings of order, Aj33, A£ 33 ~ 10 _4:_5 . 



7 



Let us finally remark that the bilinear R-parity violation leads to a distinct prediction on 
the lepton number violating decays of the lightest neutralino x i n the future colliders. The 
mass matrix of the form (J23|) enables us to determine the relation 5|£i| ^ |^| — l&l from the 
neutrino data on the mixing angles. As the parameters £j determine also the couplings of the 
R-parity violating processes; x ~^ h W t , the above mixing angle relation can be tested in 
the decay of the neutralino whose branching ratios satisfies Br(eTF) : Br(/xW) : By(tW) = 
|£i| 2 : | 1 2 : |^3| 2 01 It is intriguing to note that future colliders can provide an indirect 
test for either scenario of the axino LDM decaying through Aai or e^. 



VI. CONCLUSION 

The axino with the mass in the 1 — 10 MeV range is a good candidate for the LDM, 
which not only constitutes CDM but also explains the observed 511 keV 7-rays from the 
galactic bulge through its decay. The desired mass of the axino can be realized in certain 
supergravity models with some special arrangements, e.g., no-scale Kahler potential, or in 
gauge-mediated SUSY breaking models. The origin of relic axinos can be either the thermal 
production from the thermal bath after reheating or the non-thermal production from the 
LOSP decays. Both require a rather low reheat temperature Tr ~ 10 — 100 TeV. 

The long lifetime of the axino is a result of the R-parity violation and the suppression 
of axino interactions with ordinary particles by the PQ scale. As is well-known, the small 
violation of R-parity by bilinear terms is an attractive option for generating the neutrino 
masses and mixing. We found an interesting fact that the same small R-parity violating 
bilinear terms can explain the observed 511 keV 7-rays as well as the observed neutrino mass 
matrix consistently within the current observational bounds and the reasonable choice of 
model parameters. This connection has a virtue that the explanation of neutrino masses and 
mixing by R-parity violating bilinear terms has testable predictions in the future colliders, 
thereby provides an indirect test of decaying axino LDM. The LDM is an very attractive 
idea in that if it turns out to be true, the morphology of 511 keV gamma-rays will serve as 
a good probe of the dark matter halo density profile. The decaying LDM models require 
more curspy density profile to fit the observed morphology of 511 keV 7-rays from the 
galactic bulge than the annihilation models. We expect this leads to interesting astrophysical 
implications 31. 32l . 
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